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ABSTRACT 

Stochasticity of magnetic field lines is important for particle transport properties. Magnetic field 
lines separate faster than diffusively in turbulent plasma, which is called superdiffusion. We discovered 
that this superdiffusion is pronouncedly asymmetric, so that the separation of field lines along the 
magnetic field direction is different from the separation in the opposite direction. While the symmetry 
of the flow is broken by the so-called imbalance or cross-helicity, the difference between forward and 
backward diffusion is not directly due to imbalance, but a non-trivial consequence of both imbalance 
and non-reversibility of turbulence. The asymmetric diffusion perpendicular to the mean magnetic 
field entails a variety of new physical phenomena, such as the production of parallel particle streaming 
in the presence of perpendicular particle gradients. Such streaming and associated instabilities could 
be significant for particle transport in laboratory, space and astrophysical plasmas. 
Subject headings: MHD - turbulence 



1. INTRODUCTION 

Astrophysical plasmas feature a huge separation be- 
tween the energy containing scale and the dissipation 
scale. Such a high- Reynolds number flows a re necessarily 
turbulent (see, e.g.. lArmstrong et al.lll995j ). Conductive 
turbulent fluids generate their own magnetic fields by 
dynamo, as a result of which most of the astrophysical 
objects are magnetized to some degree, typically close to 
equipartition between kinetic and magnetic energy. One 
of the consequences of turbulence, whether in laboratory, 
astrophysical or laboratory plasmas, is the magnetic field 
line stochasticity, which plays a crucial role in most key 
physical processes, such as thermal conduction, recon- 
nection, particle transport, etc. Self-excited turbulence 
in tokamaks results in anomalous transport of particles 
perpendicular to the field and prevents reaching higher 
temperatures and densities. Recently, unusual features 
in the angu lar distribution of the arrival directions of 
cosmic ra ys (lAbbasi et a l. 2011) and spectral features for 
positrons (jAdriani et al.M2009f > revived studies of cosmic 
ray propagation in stochastic magnetic fi elds (see, e.g., 
iBeresnvak et al.ll2011t iKistler et afl l2012). Indeed, care- 
fully studying cosmic ray diffusion could help discrimi- 
nate between different models of positron excess, Fermi 
bubbles and other currently unexplained cosmic ray phe- 
nomena. 

Turbulence is now understood as a multi-scale phe- 
nomenon, byjarge owing to the pioneering paper of 
I Richardson! (|1926() . who studied turbulent diffusion and 
suggested that the diffusion coefficient depends of scale as 
D ~ I 4 / 3 , known as Richardson's law. Indeed, if two par- 
ticles are separated by a distance I, the typical separation 
speed corresponds to the typical turbulent velocity on 
scale I, which is, approximately, Sv ~ I 1 ' 3 (|Kolmogorovl 
Il941h . This suggests that separation between particles 
grow as 51 ~ t 2 ! 3 . If turbulence is uniform and character- 
ized only by the dissipation rate per unit mass, e, which 
has units of cm 2 /s 3 , it is natural that the separation be- 
tween two particles moving with the fluid, Ar, conforms 
to (Ar) 2 = goet , where go is a dimensionless number 



known as Richardson's constant. Richardson's diffusion 
has been studied extensively by experimental, theore tical 
and numerical means (see, e.g.. iSawford et af1l2008l and 
refs therein) . The turbulent diffusion of particles embed- 
ded in the MHD fluid received relatively less attention, 
however, the same type of diffusion is expected in per- 
pendicular direction due to the sam e scaling Sv ~ l 1 ^ 3 of 
strong MHD turbulence (see , e.g., iGoldreich fc Sridharl 
[19951: \Bems^v&2\2m[ [2012ah . 

A different and very interesting question is how the 
magnetic field lines separate from each other in such 
an environment. This question is crucial because well- 
magnetized plasmas are often poorly collisional, with 
the ion Larmor radius being many orders of magnitude 
smaller that the mean free path from Coulomb collisions. 
In particular, the magnetized solar wind features mean 
free paths which are comparable to the distance from the 
Sun. In galaxy clusters the Coulomb mean free path is 
around 10-100kpc. In most astrophysical environments 
there is also a high-energy component, called cosmic rays, 
for which Coulomb collisions are essentially negligible. 
Charged particles will, therefore, move along magnetic 
field lines for great distances and scatter mostly by mag- 
netic perturbations. This will result in parallel diffusion 
being much larger than the perpendicular diffusion. In 
the absence of perpendicular momentum, particle will 
move along magnetic field lines and diffuse only due to 
the magnetic field line diffusion. The motion of the bulk 
of the plasma Sv that causes ordinary diffusion could be 
neglected if the ion speed t)j is much larger than Sv □• 
For cosmic rays this condition is also very well satisfied, 
because they move along field lines with the speed com- 
parable to the speed of light c. In other words, at least for 
short timescales, the fluid is frozen from most particles' 
perspective. 

Despite being collisionless, plasmas in many circum- 

1 This is equivalent to the condition that the motions in the 
inertial range are subsonic. While some of the astrophysical tur- 
bulence feature supersonic motions on the outer scale, the inertial 
range motions are normally subsonic. 
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TABLE 1 

Simulation parameters 



Run 


Resolution 


e+/e~ 


w + /w 


e \\j e i* 


9m/ 9m 


Bl 


1536 3 


1 


1 


i 


l 


11 


512 • 1024 2 


1.19 


1.16 


1.07 


1.03 


13 


512 ■ 1024 2 


1.41 


1.37 


1.15 


1.15 


15 


1024 ■ 1536 2 


2.00 


2.36 


1.36 


1.31 


16 


1024 ■ 1536 2 


4.50 


6.70 


1.78 


1.71 



stance can be described a s fluids on scales large r 
that the ion Larmor radius (Schck o~hihin et all (2009). 
The inertial range of MHD turbulence feature strongly 
anisotropic perturbations which are much smaller in am- 
plitude that the mean magnetic field. The key compo- 
nent of this turbulence is Alfvenic mode, which is why it 
is often called Alfvenic turbulence. Due to the fact that 
the Alfven mode is driven by magnetic tension, not pres- 
sure, it is relatively unaffected by the lack of collisions. 
The presence of slow mode in suc h highly anisotropic 
turbulence neither affect dynamics (jGoldreich fc Sridharl 
[19951: IBeresnvald[20T2al) . nor influences magnetic field 
lines, as the anisotropic slow mode perturbation is mostly 
along the mean field. Therefore, the equations for the 
Alfvenic components, which are conventionally called re- 
duced MHD (RMHD), are sufficient for studying field 
lines. 

Perturbations in a strong mean magnetic field could 
be decomposed into backward and forward propagating 
eigenmodes w = v ± h/y/Airp called Elsasser variables. 
Since perturbation sources are not uniform, MHD turbu- 
lence is naturally imbalanced, i.e. the amplitudes of w + 
and w~ are not equal. This is verified by direct obser- 
vations in the solar wind, where the dom inant always 
propa gates away from the Sun (see, e.g., iWicks et al.l 
120 lit ) . Other astrophysical sources are expected to have 
strong imbalance, for example stellar winds and jets 
will emit predominantly outward-propagating compo- 
nent. Similarly, AGN jets are expected to have Alfven 
perturbations propagating away fro m the central engine, 
e.g., due to the black hole spin (jBlandford fc Znaiekl 
119771 ). The theories of imbalanced Alfvenic turbulence 
are fairly young and has been verified mostly by com- 
parison with si mulations (jBeresnvak fc Lazarianl 120091 : 
Bcresnyak 2011), although the solar wind measurements 
also show some promis e. So far the model most con - 
sistent with the data is iBeresnvak fc Lazarianl ([2008a), 
which correctly explains the ratio of anisotropies and 
the ratio of energies, given a certain ratio of energy 
fluxes. Imbalanced relativistic force-free MHD turbu- 
lence, supposedly existing in such objects as parsec- 
sc ale jets and GRB engine s has been simulated recently 
by iCho fc Lazarianl (I2013Q and seem to exhibit proper- 
ties consistent with the IBeresnvak fc Lazarianl (|2008af l 
model. Since we expect Alfven mode to survive in low- 
collisional environments such as jets, pulsar winds, etc., 
we are particularly interested in the properties of mag- 
netic field line diffusion of Alfvenic turbulence. The 
study of magnetic field diffusion is also equivalent to the 
study of charged particle diffusion in the limit of negligi- 
ble pitch angle scattering, e.g., due to zero perpendicular 
momentum. 

2. RICHARDSON-ALFVEN DIFFUSION 
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Fig. 1. — Diffusion of magnetic field lines. With initial separa- 
tions ro (0.2r;-H3.3»?) two field lines start diffusing apart at r ~ 12r; 
by Richardson- Alfven law with g m = 0.14 and transition to ordi- 
nary diffusion with coefficient 0.3Lbox a t separations around outer 
scale of turbulence. 

Assuming very strong mean field Bo pointing in the x 
direction, the equation for the magnetic field line is 



where the magnetic perturbation b = B — B is perpen- 
dicular to the mean field (Alfven mode), so that the dis- 
placement vector r will only have perpendicular compo- 
nents. Similarly, particles moving along such a field in 
one direction will only experience perpendicular diffu- 
sion, as dx 3> \dr\. Suppose, we follow magnetic field 
lines started from two points, separated by a small ini- 
tial distance ro. As the difference between B scales as 
SB[ ~ I 1 ' 3 we would expect a stochastic separation of 
the magnetic field lines in to follow the law 

((Ar) 2 )^^ 3 ^ 3 , (2) 

where e is the dissipatio n ra te per unit mass, as defined 
above, and va — B / y/4irp is the Alfven speed. This 
expression can be obtained by replacing t in Richard- 
son's formula with the time variable for the Alfven wave, 
x/va — t with t — 0. For the lack of a better term in the 
literature we will designate this Richardson- Alfven diffu- 
sion and call the dimensionless constant g m Richardson- 
Alfven constant. Although we call this diffusion by 
analogy with physical diffusion in time, this is rather 
a stochastic separation in space. Nevertheless, the term 
diffusion seems appropriate due to similarity with physi- 
cal diffusion and relevance of this problem to perpendic- 
ular diffusion of particles. 

If magnetic field lines separate for a distance much 
larger than the outer scale of turbulence L, the Sh be- 
comes truly random and independent of the separation. 
In this limit the magnetic field lines experience random 
walk, i.e. ordinary diffusion ((Ar) 2 ) ~ |a;|. This limit is 
known as field line random walk (FLRW ) and has bee n 
used to describe perpendicular diffusion in lJokipiil (|1973f) . 
Note that the random walk must be symmetric with re- 
spect to the sign of x. 

3. NUMERICAL RESULTS 

We used magnetic field snapshots obtained in simu- 
lations of Alfvenic turbulence. These simulations solved 
the reduced MHD equations with explicit dissipation and 
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Fig. 2. — The ratio of forward to backward diffusion in dat- 
acubes from simulation 15. Upper plot shows five curves from the 
same datacube with different initial separations, same as on Fig.[l] 
Lower plot shows time variability of the ratio. 

driving to achieve statistically stationary state. Fur- 
ther details behind the RMHD rationale, simulation 
setup, driv i ng, nu merical scheme, etc., can be found in 
iBeresnvakl (|2012af) . Each simulation represents station- 
ary, strong MHD turbulence with strong mean field. The 
balanced s i mulat ion Bl has been previously reported in 
IBeresn vak (201l|) and imbalanced simul a tions 11-6 has 
been reported in IBeresnvak fc Lazarianl (|2010l ) . More 
details concerning these simulations can be found in the 
above references. The parameters of the simulations are 
summarized in Table 1, with the defining feature of each 
imbalance simulation being the ratio of the dissipation 
rates for Elsasser components w . 

We were tracking the pairs of magnetic field lines 
started at random positions throughout the box and ini- 
tially separated by distance r by Eq. 1. Fig. [T] shows 
the tracking results for the Bl simulation. The tran- 
sition to Richardson's diffusion happens when particles 
are separated by around twelve Kolmogorov (dissipation) 
scales r\. We chose five initial separations, fractions of 
T). At large distances they seem to converge towards 
Richardson- Alfvcn diffusion with g m = 0.14. At suffi- 
ciently large distances the field lines started experienc- 
ing random walk, i.e. ordinary diffusion with diffusion 
coefficient of 0.3Lb ox . We typically used 4 • 10 5 field line 
pairs for statistical averaging. 

The next tracking experiment involved simulations II- 
16. These also experienced Richardson- Alfven diffusion, 
but now the diffusion speed was different depending on 
whether we track magnetic fields forward or backward 
(negative or positive dx). We plotted the ratio of for- 
ward to backward separations from the simulation 15 on 
Fig. [21 It turned out that this ratio is fairly insensitive to 
the initial separation rrj, when r$ was varied by a factor 
of around 16, with most difference being due to statistical 
error. Different snapshots of the same simulation showed 
more variation. We used this variation in time to esti- 
mate the error in the diffusion ratio. In the large x limit 
this ratio went to unity, consistent with the symmetry 
of random walk. We estimated the ratio of Richardson- 
Alfven constants by taking the maximum of the ratio 
curves, which was also somewhere around the middle of 
the inertial range in terms of perpendicular separation. 
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Fig. 3. — The ratio of forward to backward diffusion as a function 
of imbalance. 



The measurements of the ratio g^/ 1 g m are presented in 
Table 1 and Fig. El 

4. A MODEL 

Hydrodynamic turbulent diffusion forward and back- 
ward in ti me is known to b e different by a factor of a ~ 2 
(see, e.g.. lBerg et al.lr2006l ). which is due to fundamental 
non-reversibility of turbulence. In our case the diffusion 
of magnetic field lines comes from the w + component, 
which propagates against mean field and w~ component 
which propagates along mean field. The diffusion of field 
lines along the field will be "forward in time" for w~ and 
"backward in time" for w + and vice versa for the oppo- 
site direction. Assuming that the diffusion from b is 
proportional to the amplitude w ± , we can write for the 
diffusion asymmetry: 

g+_ = w~ + a m w+ 
g m a m w-+w+ , 

where a m is the time-asymmetry in MHD turbulence. 
Fig. [3] indicates that this model agrees with data reason- 
ably well, as long as a m = 2.0 -f- 2.1, which is compatible 
with hydrodynamic time-asymmetry. 

5. RELATION TO GOLDREICH-SRIDHAR ANISOTROPY 

Eq. 2 looks similar to the critical balance of 
IGoldreich &; Sridharl (|1995f l . Introducing anisotropy con- 
stant C A , critical balance can be written as £\\ = 

C A VAl 2 / 3e ~ 1/3 (jBeresnvaldl2"0T2al) . where £y is the dis- 
tance paralle to the field, i.e., analogous to x and £± 
is a perpendicular distance. This is indeed the same 
functional dependence as Eq. 2. However, if the as- 
sume that these relations are identical, this would im- 

— 1/3 

ply that Ca = gm ■ This is not satisfied, however. 
The difference is that the measurement of the diverg- 
ing magnetic fields is quasi-Lagrangian, while the mea- 
surement of structure function that lead to anisotropy 
constant is Eulerian. This difference becomes even more 
pronounced in the imbalanced case, where each of the 
w components has its own anisotropy. The analogy be- 
tween Richardson- Alfven diffusion and critical balance 
would suggest that g + (£^) 3 = g~(£^) 3 . We presented 
the ratio of parallel scales in the middle of the inertial 



4 



range in Table 1. As we see, the expression above is not 
satisfied and the ratio of magnetic diffusions can not be 
explained by the anisotropy difference. So, despite the 
similar functional form, Richardson- Alfven diffusion has 
no direct relationship to Goldreich-Sridhar anisotropy. 

6. DISCUSSION AND IMPLICATIONS FOR PARTICLE 
TRANSPORT 

Our measurement is the first demonstration of the x 3 
superdiffusion of magnetic field lines in simulations of 
MHD turbulence. Superdiffusion of fast particles in the 
solar wind has been argu ed based on observational data 
(|Perri fc Zimba rdo 200 9]) . Superdif fusion of field lines 
has been discussed in Jokipiil (|1973lh The reconnection 
model of lLazarian fc Vishniad (|1999) also uses perpen- 
dicular superdiffusion of field lin es. The superdiffusion 
of par t icles has been a r gued in |Naravan fc Medvedevl 
(2001|); lLazarian! (f2006h ; lYan~fc Lazarianl ()2008D . how- 
ever, the asymmetric superdiffusion has not been antic- 
ipated before. Our earlier measureme nt of perpendic- 
ular diffusion using MHD simulations (jBeresnvak et al.l 
120111) has been made in the large separation limit and 
reproduced FLRW, which is symmetric. The measure- 
ments of cosm ic ray propagation i n artificial random 
fields, such as iGiacinti et al.l (|201 2) can, in principle, 
reproduce superdiffusion, but since artificial fields lack 
the time-asymmetry of turbulent fields, they can not 
reproduce asymmetric diffusion. Based on the similar- 
ity between G oldreich-Sridhar anisotr o py an d Richard- 
son's diffusion INaravan fc Medvedevl (|2001|) suggested 
that magnetic field lines separate within the Goldreich- 
Sridhar cone, however, according to the section above, 
this analogy is misleading, especially in the imbalanced 
case. Time- asymmetry of turbulence, that we confirmed 



in this Letter, have conseq uences for small-scale dynamo 
as well (|Beresnvakll2012b[ ). 

One of the consequences of asymmetric perpendicu- 
lar diffusion is an induced streaming. Indeed, if we 
consider two close magnetic field tubes, one of which 
is filled with isotropically distributed particles and an- 
other empty, the asymmetric diffusion into the empty 
tube will result in an average streaming ~ (1 — g + / g~) 
of particle's velocity. In particular, for relativistic par- 
ticles, such as cosmic rays, this will result in stream- 
ing velocity of 2c(l — <? + /.<?~)/7r, which could easily ex- 
ceed the threshold for streaming instability, va as long 
as imbalance amplitude 1 — w~ /w + exceeds (3tt/2)va/c, 
which is around 10~ 4 in the WISM. Therefore, the in- 
duced str eaming will be countera cted by streaming in- 
stability (|Kulsrud fc Pearcel I1969D. Above the thresh- 
old for turbulent damping (|Farmer fc Goldreichl l2004t 
iBeresnvak fe Lazarianl l2008bf) streaming instability will 
be suppressed and, according to the estimates in the 
above paper, the weak large-scale streaming should reap- 
pear at energies 10 11 eV and strong streaming is expected 
above 3T0 13 eV, although such energies are already heav- 
ily influenced by pitch-angle scattering. The net effect of 
the streaming instability from particles with energies be- 
low 3 • 10 13 eV will be a flux of slab waves which will 
increase the rate of pitch-angle scattering for these par- 
ticles. The modeling of this effect will be subject of a 
future publication. 
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